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Abstract 

We report the first experimental observation of ac-driven phase-locked mo- 
tion of a topological soliton at a nonzero average velocity in a periodically 
modulated lossy medium. The velocity is related by a resonant condition 
to the driving frequency. The observation is made in terms of the current- 
voltage, I(V), characteristics for a fluxon trapped in an annular Josephson 
junction placed into dc magnetic field. Large zero-crossing constant-voltage 
steps, exactly corresponding to the resonantly locked soliton motion at dif- 
ferent orders of the resonance, are found on the experimental /(V) curves. A 
measured dependence of the size of the steps vs. the external magnetic field is 
in good agreement with predictions of an analytical model based on the bal- 
ance equation for the fluxon's energy. The effect has a potential application 

low-frequency voltage standard. 
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The important role played by nonlinear excitations in the form of solitons in various 
physical systems is commonly known. However, experimental observation of dynamical 
effects produced by solitons is often difficult because real systems may be far from their 
idealized mathematical models which give rise to soliton solutions. Among perturbations 
that destroy soliton effects, dissipation is the most important one. To compensate dissipative 
losses and thus make the soliton dynamics visible, one must apply an external force that 
supports (in particular, drives) a soliton0. 

Solitons of the simplest type are topological kinks, a well-known example being a mag- 
netic flux quantum (fluxon) in long Josephson junctions (LJJs)i'i. A fluxon in LJJ can be 
easily driven by bias current applied to the junction. The motion of a fluxon gives rise to dc 
voltage V across the junction, which is proportional to the fluxon's mean velocity. Varying 
the dc bias current J, one can produce a dependence V(I), which is the main dynamical 
characteristic of LJJ. An experimentally obtained I(V) curve easily allows one to identify 
the presence of one or more fluxons trapped in LJji. 

Microwave field irradiating LJJ gives rise to an ac drive acting on the fluxon. In a 
spatially homogeneous lossy system, an ac drive may only support an oscillatory motion of 
a kink, which is hard to observe in LJJ due to the absence of dc voltage. However, it was 
predictedB that an ac drive can support motion of a kink with a nonzero average velocity u in 
a system with a periodic spatial modulation. Indeed, a moving kink passes the modulation 
length (period) L during the time L/u. If this time is commensurate with the period 2tt/uj 
of the ac drive, i.e., m(L/u) = 2tt/uj with an integer m, or 

u = m[Luj /2tt), (1) 

one may expect a resonance (of order m) between the two periodicities. In other words, a 
moving kink may be phase-locked to the ac drive. This provides for permanent transfer of 
energy from the drive to the kink, making it possible to compensate dissipative losses. The 
energy balance gives rise to a minimum {threshold) value T t ^ T of the ac drive's amplitude T, 
which can compensate the losses and support the motion of the ac-driven kink. 



A more general case, when the system is simultaneously driven by the ac field and dc 
bias current /, was considered in Ref. [5[ It was predicted that the corresponding V(I) 
characteristic has steps (constant- voltage segments) at the resonant velocities (|I|). The 
motion of the fluxon under the action of pure ac drive then corresponds to zero-crossings, 
when the steps cross the axis 1 = 0. In fact, the most straightforward way to observe the 
ac-driven motion of the soliton is through zero-crossings on the V(J) characteristic. 

A formally similar feature is known in small Josephson junctions as Shapiro stepsi: ac 
drive applied to the junction gives rise to dc voltage across it (an inverse ac Josephson effect). 
However, a drastic difference of the effect sought for in this work from the Shapiro steps is 
that they are only possible at high frequencies exceeding the junctions's plasma frequency, 
while the ac-driven motion of the fluxon can be supported by the ac drive with an arbitrarily 
low frequency. This circumstance also opens way for application to the design of voltage 
standards using easily accessible sources of low-frequency radiation, which are not usable 
with the usual voltage standards based on small Josephson junctions. Moreover, because 
the single-valued dc voltage in LJJ is controlled not only by the frequency and amplitude of 
the ac drive, but also by locally applied magnetic field (see below), the system studied here 
may be useful for designing an ac voltage standard. 

The inverse Josephson effect was also studied in detail in ac-driven finite- length linear 
junctions with reflecting edges, corresponding to zero external magnetic field!. A fluxon 
turns into an antifluxon while bouncing from the edge. In that case, the system of constant- 
voltage steps is complicated, on the contrary to a very simple set of two symmetric steps in 
the circular junction, see below. This is explained by the fact that the shuttle motion of the 
polarity-reversing fluxon/antifluxon in the linear junction is different from the progressive 
motion of the fluxon, without any polarity reversal, in the circular LJJ. 

An objective of this paper is to report direct experimental observation of the ac-driven 
fluxon motion in periodically modulated LJJs. Frequently, it is assumed that the necessary 
periodic spatial modulation along the junction can be induced by periodically changing 
the thickness of the dielectric barrier separating two superconductors. In the presence of 
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the losses and drive, the modulated LJJ is described by the perturbed sine-Gordon (sG) 
equation, 



where x and t are the length along the junction and time, measured, respectively, in units 
of the Josephson length and inverse plasma frequency, e is the normalized modulation am- 
plitude, while 7 and T are the dc and ac bias current densities, both normalized to the 
junction's critical current density. 

The harmonic modulation of the local magnitude of the maximum Josephson current, 
assumed by this model, is very hard to realize in an experiment due to the exponential 
dependence of the critical current on the thickness of the dielectric barrier. A much easier 
and fully controllable way to induce a harmonic periodic modulation is to use an annular 
(ring- shaped) LJJ, to which uniform dc magnetic field is applied in its planeB. As it was 
demonstrated experimentally, fluxons can be readily trapped in annular LJJsil. In this case, 
the sG model takes the form 



where h is (renormalized) magnetic field, and R is the radius of the ring. In the case of the 
annular junction, solutions are subject to periodic boundary conditions, (p x (x+27iR) = (fi x (x) 
and <j)(x + 2ttR) = <p{x) + 2tcN, N being the number of the trapped fluxons (in this work, 
N = 1). Comparison with the experiment shows that, unlike the model (|2|), the one @ is, 
virtually, exact!. 

We assume the spatial size of the fluxon, which is ~ 1 in the present notation, to be 
much smaller than the circumference L = 2ttR of the annular junction. Large L imposes 
an upper limit on the driving frequency uj which can support the ac-driven motion: as the 
fluxon's velocity cannot exceed the maximum (Swihart) group velocity of electromagnetic 
waves in LJJ, which is 1 in our notation, Eq. ([I]) implies that uj "z l/L. 

A different type of ac drive for fluxons in circular LJJs was proposed in Ref. [j], viz., ac 
magnetic field. In this case, the terms hsin(x/R) and rsin(u;i) in Eq. (||D are replaced 
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by a single one, ft, sin (x/R) sin(u;t), which may be decomposed into two waves traveling in 
opposite directions, (l/2)[cos(x/ R — ut) —cos(x/R+ut)]. As it was shownl, either traveling 
wave may capture a fluxon, dragging it at the wave's phase velocity ±luR. A similar model 
was proposed in Ref.0, in which the fluxon is dragged by rotating magnetic field. A difference 
of our model (which corresponds to the real experiment reported below) is the separation 
between the fields that induce the spatial modulation and ac force. The separation makes 
it possible to control the dynamics in a more flexible way. 

It is straightforward to derive an equation of motion for the fluxon in the adiabatic 
approximation, following the lines of Refs. |], [5| (£= d£/dt): 

d I £ \ nh £ a £ 

cos — — 



dt vvw^ / 4Vi^7 R \A^7 

+ ^[ 7 + rsinM)] . (4) 

For further analysis, one may assume, following Refs. [7] and [H], that, in the lowest 
approximation, the fluxon is moving at a constant velocity £o= u belonging to the resonant 
spectrum ([[]), so that = ut + R5, where 5 is a phase-locking constant. Then, the first 
correction to the instantaneous fluxon's velocity, generated by the spatial modulation, can 
be easily found from Eq. (£|), 

£ 1= (nRh/Au) (l - u 2 ) sin [{u/R)t + 5} . (5) 

The approximation applies provided that the correction (||) is much smaller than the unper- 
turbed velocity it, which amounts to Rh u 2 / (1 — u 2 ). 

A key ingredient of the dynamical analysis is the energy-balance equationB. In the model 
(3) it is based on the correction (§) to the velocityB! (while in the above-mentioned models 
with ac magnetic fields00 the approximation £= u was sufficient). In the case of the 
fundamental resonance, with m = 1 in Eq. (fj), i.e., u = Ru, the energy balance yields, after 
a straightforward algebra, 

AaRuj 7rhT 
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7T,/1 - (Ruj) 2 



8Rw 2 



(Roo) 2 cos 5. (6) 



Setting |cos<5| = 1 and 7 = in Eq. (|^) gives a minimum (threshold) amplitude of 
the ac drive which can support the fluxon's motion in the absence of the dc bias current, 
T thr = (32/tt 2 /i) aR 2 oo 3 [1 - (Ru) 2 f 3/2 . For the comparison with experimental results, the 
most important consequence of Eq. (Q) is an interval of dc bias current density 7 in which 
the phase-locked ac-driven motion of the fluxon is expected. It is produced by varying cos S 
in Eq. @ between —1 and +1: 

, , AaRuj nRT [1 - (Ru) 2 } , 

vr^l - (Ru) 2 8(Ru) 2 

Note that the size of the interval strongly depends on the driving frequency, while in the 
model with the ac magnetic fielJ it does not depend on oj at all, provided that 2tcR 3> 1. 

Experiments have been performed with Nb/Al-AlO^/Nb Josephson annular junction 
with the mean diameter 2R = 95 /im and the annulus width 5 /im, applying the bias current 
i" and measuring the dc voltage V across the junction. The distribution of the bias current 
was uniform, which was concluded from measurement of the critical current I c in the state 
without trapped ffuxons at H = 0. J c was found to be about 0.9 of its value for the small 
junction. The annular L JJ had the Josephson length A j ~ 30 /im and plasma frequency 
50 GHz. Note that these parameters imply the ratio ~ 10 of the junction's length 2ttR to 
the fluxon's size, which is ~ Aj, i.e., the junction may indeed be regarded as a long one. The 
measurements were done at the temperature 4.2 K, using a shielded low- noise measurement 
setup. The ac driving current with the frequency / = u/27r between 5 and 26 GHz was 
supplied by means of a coaxial cable ending with a small antenna inductively coupled to 
the junction. The antenna was oriented coaxial to the dc bias current supplied through the 
electrodes, therefore the ac magnetic field was perpendicular to the dc magnetic field. Thus, 
the magnetic component of the ac signal induced the driving force of the same type as the 
bias current. The ac power levels mentioned below pertain to the input at the top of the 
cryostat. 

Following Ref. |^, trapping of a fluxon in the junction was achieved by cooling the sample 
below the critical temperature T c m 9.2 K for the transition of Nb into the superconductive 
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state, with a small dc bias current applied to the junction. At H = 0, the fluxon depinning 
current I^ep was found to be very small, less than 1 of the Josephson critical current I c 
measured without the trapped fluxon. As a fluxon can only be trapped by junction's local 
inhomogeneities in the absence of the magnetic field, this indicates at fairly high uniformity 
of the junction. At low values of the field H, linear increase of Jd ep with H was observed, 
which is well described by the theoretical model based on Eqs. (§) and (|j): the zero- voltage 
state exists as long as the maximum fluxon's pinning force exerted by the field-induced 
potential remains larger than the driving force induced by dc bias current, which is satisfied 
at |7| < h . Note that the fluxon depinning and re-trapping in weak external magnetic fields 
were studied experimentally and analytically in Ref. |^. 

An evidence for the progressive ac-driven motion of the fluxon is presented in Fig. |l|a. 
This I(V) characteristic was measured at H = 0.35 Oe and / = 18.1 GHz. Its salient feature 
is two large symmetric constant-voltage steps. The points where they intersect the zero- 
current axis correspond to the fluxon moving around the junction with a nonzero average 
velocity at zero dc driving force. Another remarkable feature is the absence of any step 
at the zero voltage, i.e., in the present case the fluxon cannot be trapped by the effective 
potential, even when the dc bias current is small. For comparison, in Fig. |T]b we show the 
/(V) curve measured at the same power and frequency of the drive, but with H = 0. In 
this case, a substantial zero- voltage step is seen, extending up to the current I « ±0.1 mA. 
In the absence of the ac drive, the critical current Jo is much smaller, less than 20 /iA (this 
residual Iq may be explained by small inhomogeneities of LJJ, see above). 

The conspicuous zero-voltage step in Fig. |l]b may be explained by the fact that the 
magnetic component of the ac drive creates its own modulated potential. This argument 
also helps to explain two symmetric constant- voltage steps at V ~ 37 /iV in Fig. |I|b as 
resonant steps supported by the ac-drive-induced modulation. Note, however, that the 
latter steps do not feature zero-crossing. All the data collected in the experiments show 
that the zero crossing is possible solely on the resonant steps that occur in the presence 
of dc magnetic field. In other words, the ac-driven motion of the fluxons is not possible 
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without a stationary spatially periodic potential. This inference is in no contradiction with 
numerical results of Ref. [7], where the drive itself was spatially modulated. 

Coming back to the resonant steps induced by the dc field H, which is the main subject 
of the work, we have also measured their size vs. H, see Fig. [| The result is that both edge 
values Ii and If" indicated in Fig. [I] vary nearly linearly with H, up to H pa 0.37 Oe. At 
still larger fields, the phase- locked ac-driven motion of the fluxon gets interrupted in some 
current range (the perturbation theory does not apply to so strong fields). These findings 
are in reasonable agreement with the theoretical prediction (|7|), a fit to which is shown by 
the dashed lines. This pertains to both the upward shift of the lines I t (H) (recall that / and 
H are proportional to 7 and h, respectively) and their linear change with the magnetic field. 
The residual nonzero value of I± — If at H = matches the small non-zero-crossing step in 
Fig. [l|b. It is noteworthy too that the current range of the zero-voltage state, Iq < I < Iq 
(see Fig. |]b), decreases nearly linearly with H, disappearing at H pa 0.09 Oe. 

Equation (0) also predicts a linear dependence of the step's size on the ac-drive's ampli- 
tude T. Comparison with experimental data shows an agreement in a broad power range; we 
do not display detailed results of the comparison, as they display no interesting features, the 
linear dependence of the range of existence of the inverse Josephson effect vs. the amplitude 
of the ac signal being a common feature of all the manifestations of this effect, including the 
Shapiro steps in small junctions!. 

As for the dependence of the step's size on the ac-drive's frequency at a fixed value of 
its power, it is hard to measure it, as variation of the frequency inevitably entails a change 
in the ac power coupled to the junction. Nevertheless, basic features reported in this work, 
i.e., the zero-crossing steps at finite voltages and disappearance of the zero- voltage state, 
have been observed in a broad range of the ac frequencies, starting from about 5 GHz 
and up. On the other hand, as it was mentioned above, the condition that the moving 
fluxon cannot exceed the Swihart velocity c sets an upper cutoff for the frequency that can 
support the phase-locked motion of fluxons. In our system, c corresponds to the dc voltage 
V = c$o/(2vri?) pa 80 /iV, which translates, via Eq. ([[]), into the cutoff frequency ~ 40 GHz 
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for the case of the fundamental resonance. 

All the above results pertained to the fundamental resonance, m — 1 in Eq. (^). It is also 
easy to observe zero-crossings corresponding to higher-orders resonances. This is illustrated 
in Fig. |3], showing V(I) curves with the resonant steps generated by both the fundamental 
and second-order (corresponding to m = 2) resonances. 

In conclusion, we have reported the first observation of ac-driven motion of a topological 
soliton in a periodically modulated lossy medium. The observation was made in an annular 
uniform Josephson junction placed into constant magnetic field. Experimentally measured 
data, such as the size of the constant-voltage step, are in good agreement with predictions 
of the analytical model. The effect may take place in a broad class of nonlinear systems and, 
in terms of the Josephson junctions, it may find a potential application as a low-frequency 
voltage standard. 
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FIGURES 

FIG. 1. Current-voltage characteristics for a fluxon trapped in the annular Josephson junction 
irradiated by the ac signal with the frequency 18.1 GHz and power P ac = — 8dBm. The dc magnetic 
field is (a) H = 0.35 Oe and (b) H = 0. 

FIG. 2. The critical values I Q and if 1 of the dc bias current, marked in Fig. |l|, vs. the external 
dc magnetic field, the dashed lines showing a fit to Eq. [?]. 

FIG. 3. Current-voltage characteristics for a fluxon in the annular Josephson junction at 
H = 0.40 Oe, irradiated by the ac signal at the frequency 10.0 GHz. The signal's power P ac is 
— 3.4dBm (solid line) and — 12.4dBm (dashed line). The constant-voltage steps on the two lines 
correspond, respectively, to the second-order and fundamental resonance in Eq. (111). 
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